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ANALYTIC FORMULATION OF THREE WEDGE ANGLE 


CONSTRAINTS IN TRAJECTORY STUDIES 
By Wayne 0. Laszlo 


1.0 SUMMARY AND INTRODUCTION 


For many earth-moon and interplanetary trajectory studies, it may 
be desirable that the terminal parking orbit achieved by a vehicle be 
inclined at some prespecified wedge angle (also called angle of inclin- 
ation) measured with respect to seme reference coordinate system 

. ■+• -*■ . , 

(i,j,k unit vectors). The occupied focus of the terminal parking orbit 

would be the center of the attracting planet. In the ehrth-moon 
trajectory optimization studies presently being undertaken by the 
Mathematical Physics Branch, Mission Planning and Analysis Division, 
a wedge angle constraint is being used in conjunction with a terminal 
circular lunar packing orbit. 

The purpose of this document is to derive an analytic formulation 
of three wedge angle constraints for a vehicle’s terminal parking orbit 
from a given state vector on that orbit. In the analysis presented, 
it is determined that the use of a particular wedge angle constraint 
for a vehicle's terminal parking orbit depends on the freedom desired 
in both the vehicle's direction of motion (posigrade or retrograde) 
on that orbit and in the orientation of the orbital plane. The three 
wedge angle constraints and their respective interpretations (including 
geometrical interpretations) are as follows. 

First, when cos W - cos W = 0 is satisfied, the vehicle's direction 
of motion corresponding to W will be fixed, and the orientation of 
the orbital plane will be left free to the extent that H rotates to 
form a one-nappe cone of angle of revolution W about the k axis. 
Lecond, when sin W - sin W = 0 is satisfied, the vehicle's direction 
of motion corresponding to W will be left free, and the orientation 
of the orbital plane will be free to the extent that H rotates to form 
a two-nappe cone of angle of revolution W about the k axis. This 
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constraint allows the maximum freedom in the vehicle's terminel parking 

orbit of the three constraints. Third, when H - (h) =0 is satisfied, 

the vehicle's direction of motion corresponding to W will be fixed, and 
the orientation of the orbit plane will be fixed to the extent that there 
is no cone of revolution. This constraint allows the minimum freedom in 
the vehicle's terminal parking orbit of the three constraints. 


cos W 
cos W 
e 
H 


|H|, H 
(H) 





M 


| M | , M 
N 


2.0 SYMBOLS 

cosine of the wedge angle 

cosine of the fixed wedge angle W 

eccentricity of the terminal parking orbit 

angular momentum vector associated with an orbit. 

The components in the coordinate system 

are H , H , H 
x' y * z 

magnitude of H (or /H z +H z + H z ) 

x y z 

fixed angular momentum vector that is to be 
achieved by sati sfac tion of the wedge angle 

constraint H - (H) = 0 

mutually orthogonal unit vectors in the reference 
coordinate system 

mutually orthogonal unit vectors in an alternate 
coordinate system 

“► “► "► 

vector given by M = N * k. The components in the 

/“► -► i 

(i,j,k) coordinate system are M , M , M 

x y Zi 

magnitude of M (or /m^ z + M^ 2 + M^ 2 ) 

vector of ascending node given by N = k * H. 

The components in the (i,j,kj coordinate system 

are N , N , N 
x’ y z 
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[ N | , N magnitude of N (or /N ^ + N 7 + N 

x y z 

R(x,y,z), V(u,v,w) state vector of the vehicle on its terminal parking 

orbit. The position components are x,y,z, and 
the velocity components are u,v,w. Components 

are in the (i,J,k) coordinate system. 


sin W 
sin W 
V 

W 

W 


sine of the wedge angle 

sine of the fixed wedge angle W 

angle determined by the dot product of M and H , 
radians 

wedge angle for any parking orbit determined by the 
dot product of k and H , radians 

fixed wedge angle that is to be achieved by satis- 
faction of the wedge angle constraint , radians 


x 

a 


ft 


denotes cross product between two vectors 

fixed angle of inclination of the ( i * , j * ) plane 
with respect to the (T,j) plane, radians 

fixed longitude of ascending node (or angle between 
i and i’)» radians 

denotes dot product between two vectors 


3.0 ANALYSIS 


The wedge angle W for any parking orbit will be defined as the 

angle determined by the dot (•) product of k and H, where H is the 
angular momentum vector associated with the orbit. The wedge angle is 

also given by the solid angle between the reference (i,j) plane and the 
parking orbit plane. By the geometry of figure 1, it can be concluded 

that 0 < W < tt. Furthermore, for 0 < W < 77 , H >0 and the orbit 

will be posigrade; while for 77 < W <_ tt , H =0 and the orbit will be 

u Z 

TT 

retrograde. If W = — , then H =0 and the orbit may be called a polar 

u Z 

orbit . 
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To derive an analytic formulation of one wedge angle constraint for 
a vehicle's terminal parking orbit, it is assumed that R(x,y,z), 
V(u,v,v) represents a state vector of the vehicle for that orbit. Let 
W (a fixed angle in radians) be the wedge angle that is to be achieved. 

Hj *► - . 

The vector H ~ R x V (* indicates cross product) is computed. The 
determinant given by equation (1) defines H. 


H = 



3 

y 

V 



lll + HJ+Hk 
x z 


(1) 


Thus,, the components cf H can be written according to equations (2), 
( 3 ) , and ( 4 ) . 


H = yw - zv 
x J 

K = zu - xw 

y 

H = xv - yu 
z 


(2) 

(3) 

(M 


If H = 0 (i.e.,H = H = H = 0 ) , then R is parel] el to V, which 

x y z 

implies rectilinear motion and impact at the center of the planet . 
Although such an event is very unlikely, the wedge angle constraint will 

diverge if H = 0, and a numerical value of the constraint would have to 
be specified. 


If H # 0, cos W can be computed from equation (5). 

— = cos W (5) 

1*1 |h| 

From the latter equation, note that W is the smallest angle between 
k and H and that 

H 

cos W = jp (6) 

where H = |H| . Then, given W and therefore cos W, the wedge angle 
constraint may be written as 


cos W - cos W = 0 


( 7 ) 
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In terms of (P,V), 


equation (7) may be written as follows. 


xv - yu 

/H y + H ^ + H ^ 
x y z 


cop W = 0 


( 8 ) 


where H , H , H 
x y z 

equations (2), (3), and (U). __ _ 

[given by equations (7) or ( 8 ) J is precisely equivalent to the constraint 


are given in terms of x, y, z, u, v, w by 

Note that because 0 < W < tt , this constraint 


W - W = 0. To give a geometrical interpretation of the wedge ingle 
constraint [eq. (7)], recall that the wedge angle W for an orbit carries 
with it the vehicle's direction ol motion (posigrade or retrograde) in 

the sense that 0 <_ W < ^- implies posigrade motion while ~ < W <_ tt 

implies retrograde motion. Then, because equation (7) is equivalent to 

the constraint W - W = 0 (W an input constant), it follows that when 

cos W - cos W = 0 the vehicle will nave achieved a terminal parking 


orbit with the fixed wedge angle W as well as with the fixed direction 

t i 

of motion corresponding to W. However, while equation (7) does constrain 

-*• -► -► 
the angle between k and H, namely W, the direction of H has some 

freedom of rotation. In other words, the analytic formulation given 

ny equation (7) yields some freedom of orientation of the terminal 

parking orbit plant. In fact, the geometrical interpretation of 

equation (7) (fig. 3) is that H rotates to form a one-nappe cone of 

angle of revolution W about the k axis. The freedom of orientation 
of the orbital plane is also shown in figure 3 in that orbits 1 and 2 
each have permissible planar orientations which satisfy equation (7). 

In trajectory optimization studies in which equation ( Y ) is used, the 
optimal (i.e., minimum fuel) orientation of the terminal parking orbit 
plane will be attained. 


Suppose it is desired to derive an analytic formulation of a wedge 
angle constraint which when satisfied will leave free the vehicle's 

direction of motion corresponding to W (posigrade or retrograde). 

Such a constraint would result in the same freedom of orientation of 
the orbital plane as equation (7). The constraint that will permit 
satisfaction of the latter two functions is given by equation (9). 


sin W - sin W = 0 


(9) 
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For 0 <_ W <_ it, it follows that if sin W * sin W (W an input constant), 

then W=WorW= 7 r-W. Therefore, if 0 <_ W < ^7 (posigrade motion), 

then conceivably W * it - W (retrograde motion). It may then be said 
that when equation (9) is satisfied, the vehicle's direction of motion 

corresponding to W is free, and there is the same freedom of orientation 
of the orbital plane as in the previous formulation. 


To derive the analytic formulation of the wedge angle constraint 

[eq. (9)], from R(x,y,z), V(u,v,w), H is computed by equation (l) 
with components given by equations (2), (3), and (H). Then, if 

H f 0, N * k x II, the vector of ascending node is computed. The vector 

N is given by the determinant 



3 

0 

H 

y 




+ H J + N k 

y z 


do) 


Thus, the components of N can be written according to equations (ll), 
( 12 ), and (13). 


N = -H 

( 11 ) 

x y 

N = H 

( 12 ) 

y x 

N = 0 

rj 

(13) 


If N = 0 (i.e., -H = H = 0 , or k is parallel to 

y x 

is obtained with either a C wedge angle (inclination) 
tt wedge angle (inclination) if H <0. However, if 
either H j 0 or H ^ 0 or both), then M = N x k 

y 

The vector M is given by the determinant 


H), then an orbit 
if H > 0 or with a 

- I 

N i 0 (i.e. , 

is computed. 


M = det 


/I 

("■ 

v 0 


J 

N 


0 = M 1 + M J + Mk 

1 x y° z 


(1U) 


0 


1 - 
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f 


f 


Thus, the components of M can be written according to equations (15) » 
(16), and (17). 


M 


M 


N 


y 

N 


M - 0 
z 


H 


x 

H 


(15) 

(16) 
(17) 


O 


Next, it can be shown that M is a vector in both the (i,J) plane 

•+ -t 

and in the plane determined by the noncollinear vectors k and K. To 

*♦ 

verify this reeult, recall that K must be perpendicular to N and k 

from vector analysis because M = N x k. Further, because N = k * H, 

it car. be concluded that all vectors in the plane determined by vectors 

k and H are perpendicular to N, and any vector perpendicular to N 

must lie in this plane. Therefore, M lies in the plane determined by 

vectors k and H. With the latter results in mind, the angle V will 

be the angle determined by the dot (•) product of M and H, that is, 

the smallest angle between M and H. As seen in figure 2, if 0 < W <_ 

then 


W ♦ V = 2 


(16) 


Further, if — < W < it , then 


V - W - £ 


(19) 


For either case [eq. (18) or eq. (19)] , the cosine of V can be written 
as follows . 


cos V = 


M 


H 


— ► *T 

m| |h 


= sin W 


wnere 


c 


and 


Ml = M = 

x y z 


H | = H = /H^ + H^+H 2 
1 x y z 


(20) 


o 


ro|^ 
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The rewriting of equation (20) results in equation (21) 

MH ♦ M H + MH 

. fin W * x x 5-2- (21) 

and the use of equations (15 )» (lo), and (if) results in equation (22). 

H 1 ♦ H 2 

sin W * x (22) 

✓H 2 ♦ H 2 + 

x y x y 7. 

Therefore, in terms of (R,V), equation (9) may be written as 

K 2 + H 2 

■ - ■ * ■ y - ■■■ - sin W = 0 (23) 

✓H ' l * H 2 /H * ♦ H ' l H * 
x y x y z 

where H , H , H are given in terms of x, y, z, u, v, w by 
x y z 

equations (2), (3), and (U). 

The geometrical interpretation of equation (9) (fig* *0 is that H 

rotates to form a two-nappe cone of angle of revolution W about the 

k axis. The freedom of orientation of the orbital plane is shown in 
figure U in that orbits 1 and 2 each have permissible planar orientations 
which satisfy equation (9). 

Finally, suppose it is desired to derive an analytic formulation of 
a wedge angle constraint which when satisfied will fix the orientation 
of the orbital plane and will fix the vehicle's direction of motion 

corresponding to W. The constraint that will result in satisfaction of 
the latter two functions is given by 

H - (IT = 0 (2U) 

where H is given in terms of x, y, z, u, v, w by equations (2), (3), 
and (U). 

The geometrical interpretation of equation ( 2h ) is given in 
figure 5* As seen by the figure, there is no cone of revolution to 
correspond to equation (2U), and orbit 1 has the only permissible planar 
orientation that satisfies equation (2U). 
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U.O CONCLUDING REMARKS 

Three comments are in order. First, in the event that W = 0, 7 -, 
or ir , two of the three constraints are equivalent. For W = 0, 

-► ,-+■ . — - IT 

H - (H) = 0, and cos W - cos W = 0 are equivalent; for W = — , 

cos W - cos W = 0 and sin W - sin W = 0 are equivalent; and for 

W = tt , H - (H) = 0 and cos W - cos W = 0 are again equivalent. Second, 
to insure that the vehicle's terminal parking orbit is either elliptical 
or circular, a constraint of the form 0 <_ e < 1 (e denotes eccentricity 
of the terminal orbit) may be used in conjunction with the desired wedge 
angle constraint. Third, note that the wedge angle W for any parking 

orbit is measured with respect to some reference (i,j,k unit vectors; 
coordinate system. However, suppose it is desired that the vehicle 

achieve a parking orbit with wedge angle W measured with respect to 

some alternate coordinate system (T',j^,k' unit vectors). If the vehicle 

is assumed to be located at some state vector R(x,y,z)', V(u,v,w) in 

* -> -> “► . 

space, it will be necessary to rotate the (i,j,k) system into the 

.■f -► — ► . "► 

(i',j',k') system and to express the state vector R ,V in terras of the 
. “► -> •> . 

alternate (i',j',k') system. The latter process could consist of two 
rotations, the first through the angle 51 from 1 to T' (51 called 
the longitude of ascending node ) and the second through the angle a 
from the (i,j) reference plane to the (T',j') reference plane [a called 
the angle of inclination of the (i',j') plane with respect to the 
( i , j ) plane ] . 


v 
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IT 



it 


0 VEHICLE'S TERMINAL PARKING ORBIT 

SV TYPICAL STATE VECTOR ON ORBIT 0 
ft THE VECTOR OF ASCENDING NODE 

it THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH POSIGRADE 
MOTION ON ORBIT 0 

W THE WEDGE ANGLE ASSOCIATED WITH POSIGRADE MOTION ON 
ORBIT 0 

H* THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH RETROGRADE 
^ MOTION ON ORBIT 0 

(ir - W) THE WEDGE ANGLE ASSOCIATED WITH RETROGRADE MOTION ON 

ORBIT 0 


Figure 1.- Geometrical interpretation of the wedge angle in reference coordinate system. 
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0 VEHICLE'S TERMINAL PARKING ORBIT 

SV TYPICAL STATE VECTOR ON ORBIT 0 

N* THE VECTOR OF ASCENDING NODE GIVEN ByIT^xIT 
N? THE VECTOR GIVEN BY ^InTxT 

IT THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH POSIGRADE 
MOTION ON ORBIT 0 

W THE WEDGE ANGLE ASSOCIATED WITH POSIGRADE MOTION ON 
ORBIT 0 

V THE ANGLE DETERMINED BY THE DOT (•) PRODUCT OF M*AND H* 

P PLANE CONTAINING VECTORSXTL ANDM* 


Figure 2.- Vector geometry. 


12 


) 




ORBIT WITH ASSOCIATED ANGULAR 
MOMENTUM VECTOR^ 

ORBIT WITH ASSOCIATED ANGULAR 
MOMENTUM VECTOR"^ 



H, THE ANGULAR MOMENTUM VECTOR 
1 ASSOCIATED WITH POSIGRADE . 
MOTION ON ORBIT 1 

H 9 THF ANGULAR MOMENTUM VECTOR 
* ASSOCIATED WITH POSIGRADE 
MOTION ON ORBIT 2 


W THE FIXED WEDGE ANGLE THAT 

IS TO BE ACHIEVED BY SATISFYING 
THE WEDGE ANGLE CONSTRAINT 
COS W - COS W= 0 


Figure 3.- Geometrical interpretation of the wedge angle constraint cos W - cos W = 0. 
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1 ORBIT WITH ASSOCIATED ANGULAR MOMENTUM VECTOR EITHER H* OR H^ 

2 ORBIT WITH ASSOCIATED ANGULAR MOMENTUM VECTOR EITHER OR H^ 

H^ THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH POSIGRADE MOTION ON ORBIT 1 

THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH RETROGRADE MOTION ON ORBIT 1 

H^ THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH POSIGRADE MOTION ON ORBIT 2 

H* THE ANGULAR MOMENTUM VECTOR ASSOCIATED WITH RETROGRADE MOTION ON ORBIT 2 

W THE FIXED WEDGE ANGLE THAT IS TO BE ACHIEVED BY SATISFYING THE WEDGE ANGLE 

CONSTRAINT SIN W- SINW= 0 



Figure A. - Geometrical interpretation of the wedge angle constraint sin W - sin W = 0 
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ORBIT WITH ASSOCIATED ANGULAR MOMENTUM VECTOR 

THE FIXED ANGULAR MOMENTUM VECTOR THAT IS TO BE 
ACHIEVED BY SATISFYING THE WEDGE ANGLE CONSTRAINT 

H - |gj = 0 

THE FIXED WEDGE ANGLE THAT IS TO BE ACHIEVED BY 
SATISFYING THE WEDGE ANGLE CONSTRAINT 

H - = 0 



Figure 5.- Geometrical interpretation of the wedge angle constraint H 
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